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George Kotsopoulos 
1906-373 Front St. West, Toronto, ON M5V 3R7 

(kq ; g.kotsopoulos@utoronto.ca 

7— I 
O 

(N ; and 

a 

Charles C. Dyer 

in; 

Dept. of Physical and Environmental Sciences, University of Toronto Scarborough, and 
Dept. of Astronomy and Astrophysics, University of Toronto, Toronto, Canada 

u 

tJD dyer@astro .utoronto . ca 

January 24, 2012 

m 
in 

o 

(N 



X 



Received ; accepted 



-2 - 
Abstract 

Plane symmetric perturbations are applied to an axially symmetric Kasner 
spacetime which leads to no momentum flow orthogonal to the planes of symme- 
try. This flow appears laminar and the structure can be interpreted as a domain 
wall. We further extend consideration to the class of Bianchi Type I spacetimes 
and obtain corresponding results. 

Subject headings: general relativity, Kasner metric, Bianchi Type I, domain wall, 
perturbation 
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1. Finding New Solutions 

There are three principal exact solutions to the Einstein Field Equations (EFE) 
that are most relevant for the description of astrophysical phenomena. They are the 
Schwarzschild (internal and external), Kerr and Friedman-LeMaitre- Robertson- Walker 
(FLRW) models. These solutions are all highly idealized and involve the introduction of 
simplifying assumptions such as symmetry. 

The technique we use which has led to a solution of t he EFE involves the perturbation 



of an already symmetric solution ( IWilson and Dyer 



20071 ). Whereas standard perturbation 



techniques involve specifying the forms of the energy-momentum tensor to determine the 
form of the metric, we begin by defining a new metric with an applied perturbation: 

9ab = 9ab + Kb (1) 

where g a b is a known symmetric metric and h a b is the applied perturbation. We then use 
this new metric to find the components of the Einstein tensor, G a b, in terms of the metric 
components. Next, we invert the EFE, where now T a b = G a b/n, to ascertain whether or not 
components of the energy-momentum tensor, T a b, exist to satisfy the perturbation. If such 
an energy-momentum tensor is physically acceptable, we then try to limit the behaviour 
of our perturbation by imposing conditions to model particular astronomical phenomenon. 

te structure, we could choose our density 



2007J). This makes the 



For example, if we were modelling a galaxy-lir 
to scale as p ~ r~ 2 from the galactic centre (IWilson and Dyerl 
technique flexible and relevant in that we can model a physically viable structure. With 
these real physical constraints applied, we have a new metric. We now re-compute the 
energy-momentum components of this new metric and investigate any implications that 
may arise from our new exact solution to the EFE. 
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The Base Solution 



The metric we choose to perturb is the Kasner spacetime (lKasnerlll92l[ ). This metric 
describes an anisotropic universe that is a vacuum solution to the EFE classified as a 
Bianchi Type I universe. It is characterized by a spacetime that is anisotropically expanding 
(or contracting) in two directions while contracting (or expanding) in the third with 
space-like slices that are spatially flat and with a singularity occurring at t = 0. Hence, 
the Kasner spacetime is of special interest in cosmology since the standard cosmological 



solutions to the EFE near the cosmological singularity, suc h as the aforementioned 



model, can be described as a succession of Kasner epochs ( ILifshitz and Khalatnikov 



TRW 



19631 ) 



Furthermore, by breaking a homogeneity along one directio n, Bianchi Typ e solutions can 



easily be generalized to spacetimes with G2 isotropy groups ( jHarvey 



19901 ). 



The Kasner spacetime in Cartesian coordinates with metric signature 
takes the form (in spacetime dimension D): 

D-1 

ds 2 = dt 2 -j2 t2P >[ dx j} 2 

where P~ are the Kasner exponents which satisfy: 



j ■>•••■> 



(2) 



D-1 



D-1 



= 1 and J2 P I 



1 



(3) 

3=1 3=1 

The first condition in (3) describes a plane whereas the second condition describes a sphere 
of dimension D — 1. Thus, the Kasner exponents lie on a sphere of dimension D — 2. For 
D = 4, at each t = constant hypersurface, there exists a flat 3-dimensional space, whose 
worldlines of constant x, y and z are time-like geodesies along which galaxies, or other 
test particles, can be imagined to move. Since the solution represents an anisotropically 
expanding (or contracting) universe, a volume element dV increases (or decreases) in time 
as A/]gj<i 3 3; = td 3 x, where g = det\g a b\. 
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The Kasner spacetime in Cartesian coordinates in a 4 — D spacetime takes the form: 

ds 2 = dt 2 _ t 2P ldx 2 _ t 2P 2dy 2 _ t 2P 3dz 2 ^ 

with the same restrictions imposed by the Kasner exponents as given in (3). 

The form of the Kasner spacetime that we consider is the axisymmetric case with G2 
isometry that occurs when two of the exponents are equal, which by (3), requires that 
Pi — P 2 — 2/3 and P 3 = —1/3. With this restriction, the metric in cylindrical coordinates 
becomes: 

rfs 2 = df 2 _ t 4/3 (rfr 2 + r 2 rf02) _ t -2/3 dz 2 (5) 

3. The Perturbations 

The set of perturbations we implement are only on the ^-component of the Kasner 
spacetime represented in cylindrical coordinates (5). Clearly, the results we obtain in 
cylindrical coordinates will be similar to those in Cartesian coordinates since r 2 = x 2 + y 2 . 
For simplicity, we define the Kasner spacetime in cylindrical coordinates as: 

ds 2 = dt 2 - t 4 / 3 (dr 2 + r 2 d<p 2 ) - g zz dz 2 (6) 

where g zz = g zz + h a b is our perturbed metric. We consider the three Cases shown in Table 
1. 



Cases 


9zz = 9zz + Kb 


I 


[r 2 / 3 + h] 


II 


[r 1 /- 3 + H] 2 


III 


[t- 2 / 3 + Z(z)T(t)\ 



Table 1: Implemented perturbations. 
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In Cases I and II, we first investigate H as a different iable function of z. We then investigate 
Cases I and II again, but where H is now made a differentiable function of t, r and z. In 
Case III, we investigate the result of a perturbation with a product separable differentiable 
solution. 



To aid in these computations, we use the REDU CE Computer Algebra System (IHearn 



20091 ) with the REDTEN Tensor Analysis Package ( IHarper and Dyer 



19941 ) 



4. Results of the Perturbations 

Both Cases I and II produce a G a b with non-zero components. When the perturbation 
is H(z), the result produces non-zero terms for G a , G rr and G^. All other terms, including 
G zz , are zero. When the perturbation is H(t,r,z), a similar result is obtained for all 
cases, but with the addition of cross-components, G a b for a ^ b. These cross terms can 
be made zero if we choose the perturbation to contain only first-order terms in H or by 
making OH /Or = dH/d<f) = 0. Also, in cylindrical coordinates for all the cases reviewed, 
G rr = r 2 G </,</,, which is expected. The product separable perturbation of Case III reveals 
the same results as Case I and II; we arrive at non-zero terms for Gu, G rr and G^ and a 
G zz = 0. Thus all perturbations that involved only the g zz term of the Kasner metric reveal 
that we will always obtain G zz = when considering axisymmetric or plane symmetric 
Kasner spacetimes, when PI = PI. 

Relating these results to the energy-momentum tensor T a b, both the G zz and G za = G az 
(for a = {0, 1,2}) components of the Einstein tensor being zero reveals that there is no 
momentum-flux across the z = constant surface regardless of the type of perturbation 
applied. Each of the applied perturbations resulted in the metric collapsing (or expanding) 
in the z-direction while expanding (or collapsing) in the x- and y-directions, as is expected 
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of the Kasner metric. But, as a result of the G zz term being zero, there is no interaction 
between the stratified layers of the matter above or below the xy-plane. This is suggestive 
of laminar flow. 

Inverting the EFE to get T a b = G a b/n, leads to T ao being diagonal but with no T zz 
component. The energy-momentum tensor for an infinite, static plane-symmetric domain 
wall as suggested by Campanelli et al is: 



Tab = S(z)diag(p, -p, -p, 0) 



(7) 



whe re p is the energy-dens ity and p is the pressure. This energy-momentum ten- 



sor (ICampanelli et al 



wall ((Kibble 



20031 ) corresponds to an infinite, static plane-symmetric domain 
19761 ) lying in the xy-plane. Domain walls correspond to a particular class 
of topological defects whereby the energy-density is trapped, and from a cosmological 



point of view, these over-dense regions cou 



d lead to structure fo rmation (iBrandenberger 



19971 ). Furthermore, it has been suggested (IFriedland et al.l 120031 ) that the Universe may be 



dominated by a network of domain walls and these domain wall structures could represent 
an alternative view of dark energy theories. 



5. Energy Conditions 

In order to rule out any non-physical solutions to the EFE, energy conditions are 
applied to the state of matter content for gravitational and non-gravitational fields. These 
energy conditions consist of the Weak, Null, Strong and Dominant energy conditions, which 
are coordinate-invariant constraints on the energy-momentum tensor. 

For the cases in which our Kasner spacetime was perturbed with H(z), we apply the 
least stringent of these conditions, the Null Energy Condition (NEC). This condition states 
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that for all future-pointing null vectors, k a : 

T ab k a k b ^ 



(8) 



This restriction implies that p + p ^ 0, whereby the energy-density may be negative as long 
as there is a compensating pressure. If this condition is violated in any of the perturbed 
Kasner spacetimes, it would then indicate that our solutions are unstable. 



We choose the null vector field k a , in the Kasner spacetime to be, 

k a = 



1 - n 2 g z 



9tt 



-,o,ox 



(9) 



where neM.. This class of null vectors lies in the T zz plane of interest. A summary of the 
energy-density restrictions is shown in Table 2, where T = dT/dt. For all Cases I to III, 
T ab k a k b is positive definite for positive H. 



Cases 


Qzz — Qzz ~\~ h a \) 


T ab k a k b 


I 


[t- 2 / 3 + H(z)\ 


4fV 3 (t 1 / 3 n 2 +t)if 
9t 4/3( 6t 2/3# +1 ) 


II 


[t- 1 / 3 + H(z)} 2 


A(t l l Z n 2 +t)H 
9i 4 / 3 (t 1 /3 i j + i) 


III 


r 2 / 3 + H{z)T{t) 


2t 1 / 3 (3t 1 / s fn 2 t+3Tt 2 +2t 1 /*Tn 2 +2Tt)H 
9t 4/3( 6t 2/3# T+1 ) 



Table 2: Applied NEC results 



6. Bianchi Type I Metrics 



As mentioned in Section [21 the Kasner spacetime is a special class of Bianchi Type I 
spacetime being a homogeneous and anisotropic vacuum solution to the EFE. Therefore, 
it is appropriate to apply the perturbations to the more generalized Bianchi Type I to 
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investigate its d ynamics and co mpare the results to the perturbed Kasner models. The 



Bianchi Type I (jStephani 



20031 ) metric has the form: 



ds 2 = c 2 dt 2 — g\\dx l — g-ndy 2, — g-^dz^ 



(10) 



where, 

g aa = (-g) 1/3 [ct/(Mct + A)] 2P -- 2 ^ (11) 

such that a = 1,2,3, with no sum on a. The quantity M = K[ic 2 ^/—g is a constant and 
yj— g = 3ct(Mct + A)/4, where A is an integration constant. Using these relationships and 
an appropriate constant rescaling of coordinates, we can re-write (TTU]) as: 

ds z = dt 2 - t^dx 2 - t^dy 2 - r 2 / 3 (l + ttfdz 2 (12) 

where e = M/A. As before, since we wish to compare this result to the axisymmetric 
perturbed Kasner spacetime, we impose Pi = P2 = 2/3 and P3 = —1/3. 

There are two limiting cases of the Bianchi Type I metric in (12). If we choose et <C 1, 
the initial singularity is approached at early times and we regain the original Kasner 
spacetime. For late times, if we choose et ^> 1, the metric isotropizes to the well-known 
Einstein-de Sitter dust metric, which is a sub-set of the FLRW metric. Indeed, the Kasner 
solution is a past asymptotic state as mentioned in Section |2j 



7. Perturbations on Bianchi Type I 



The perturbations applied to the Bianchi Type I metric are summarized in Table 3. 



- 10 - 



Cases 


9zz 9zz hab 


I 


t- 2 / 3 [{l + et) 2 + H(z)\ 


II 


t-^[(l + et) 2 + H(z)T(t)] 



Table 3: Implemented Perturbations on Bianchi Type I 

Investigating G a b for the perturbations of Cases I and II, we arrive at the same conclusion 
as was found for our perturbed Kasner spacetimes; we obtain G zz = for the axisymmetric 
Bianchi Type I spacetimes. 

8. Symmetries 

The nonlinear nature of the EFE makes it difficult to find exact solutions. All of the 
known solutions have admitted simplifying symmetries in order to attain a solution. The 
technique we apply involves the perturbation of an already highly symmetric spacetime in 
hopes of producing a new spacetime solution. We now wish to determine to what extent 
any symmetries that were inherent in the original spacetime metric still retain (or break) 
any symmetries. 

To investigate symmetries, we will consider conformal Killing vector fields, £ a , that 
satisfy the conformal Killing equations: 

£o||6 + &||a = <t>9ab (13) 

where = \t°\\ c an d II denotes covariant differentiation. Utilizing the REDUCE/REDTEN 
computer algebra system, we computed the conformal Killing equations as 10 symmetric 
rank-2 tensors for each of the perturbation cases investigated. A Killing vector from the 
original metric was then calculated along with its covariant derivative and substituted in 
the conformal Killing equation (13). For the perturbations applied to the two Cases, no 
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such Killing vectors or conformal Killing vectors were found to remain; the original Killing 
vector is no longer a Killing vector. 

9. Conclusion 

Starting from two known, highly symmetric, solutions to the Einstein Field Equations 
we have applied plane symmetric perturbations and have shown that the resulting perturbed 
spacetimes exhibit the existence of structure that can be interpreted as a domain wall. It 
was demonstrated that these solutions do not violate the Null Energy Condition and thus 
would permit us to apply an appropriate energy-momentum tensor in future investigations. 
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